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1.0 INTRODUCTION 


The concept of magnetic bearings basically consists of using magnetic forces 
to maintain controlled separation between bearing surfaces. The magnetic 
bearings offer a number of advantages over conventional bearings: 

• They require no lubricants and, therefore, can operate in vacuum 
and be free from contamination; 

• They have low mechanical friction and therefore energy can be saved; 

• The start-stop rub is absent; 

• The bearing stiffness can be controlled by electrical means, thus 
helping to step through critical speeds. 

The magnetic bearings can employ either closed- or open-type magnetic circuits. 
An example of the latter case is two magnets in repulsion. In this case, the 
repelling force increases as the two magnets are brought closer, indicating 
stability. However, at the same time, it is difficult to keep the magnets 
aligned because any small transverse displacement creates a transverse force 
which tends to push the magnets further apart. The system is thus unstable in 
the transverse direction. The stable and unstable directions are also referred 
to as passive and active directions, respectively. The closed-circuit system, 
exemplified by two attracting magnets, is active in the normal direction and 
passive in the transverse direction. 

Magnetic bearings are thus not inherently stable in all the three directions, 
and this has been shown to be true mathematically (Earnshaw's Theorem, [1]*). 
For successful operation, therefore, the active directions have to be servoed. 
Using the method of electrical servoing, many active magnetic suspensions have 
been built [2, 3]. However, maintaining passive forces as much as possible 
decreases the need to depend on electrical servoing and thereby reduces the 
electrical power requirements and propensity for instability. 

Techniques to predict magnetic forces in both open- and closed-circuit sys¬ 
tems have not yet been fully developed. However, over the past several years 


*Numbers in brackets indicate references in Section 8.0 of this report. 









a considerable amount of research on magnetic suspensions has been performed 
at Mechanical Technology Incorporated (MTI). The continued sponsorship by 
the Office of Naval Research (ONR) has resulted in the development of a vali¬ 
dated theory for predicting forces and fluxes in a closed-circuit system 
employing high-permeability ferromagnetic materials. This work has been de¬ 
scribed in an earlier report [4]. The present report contains a description 
of all the research performed for ONR on magnetic bearings subsequent to that 
reported in Reference 4. 

The research described in this report contains the results of studies performed 
on both closed- and open-circuit magnetic suspensions. The closed-circuit 
system was the same as that investigated earlier [4]. The studies performed 
on closed-circuit suspension consisted of the production of analytical curves 
for the prediction of forces and fluxes and further experimental verification 
of the theory when permanent magnets are used to energize the magnetic circuit. 
The work on open-circuit suspensions consisted of the development of a theory 
together with experimental verification for predicting forces in a quasi¬ 
static environment. These are described in more detail below. 

The theory for predicting forces in closed-circuit systems [4] was used to 
produce curves for nondimensional values of normal and transverse forces. The 
variables were gap thickness and eccentricity, and the parameter considered 
was related to tooth geometry. Results have been reported for both straight- 
and circular-tooth configurations. 

Dynamic measurements were performed on the same closed magnetic circuit used 
previously in order to determine any possible effects of transverse vibration. 
Although the measurements indicated some degradation of the transverse forces 
at moderately high frequencies of vibration, the phase-angle measurement be¬ 
tween force and displacement indicated an extremely small value of damping in 
the system. 

Towards the study of open-circuit suspensions, a theory was developed to pre¬ 
dict both normal and transverse forces between two repelling permanent mag¬ 
nets, together with its experimental verification. For purposes of analysis, 
the magnets were mathematically modeled by current densities, and forces were 










calculated using the Ampere's Law. A modified version of the test apparatus 
used for dynamic measurements was utilized for experimental measurements. 



II 

Also, the theory that was developed for closed-circuit systems employing high- 
permeability materials was verified for the case when permanent magnets are 
used as a battery to energize the magnetic circuit. The rare-earth-cobalt 
permanent magnets were used in the experimental measurements. 














2.0 ANALYTICAL RESULTS FOR CLOSED MAGNETIC CIRCUIT WITH HIGH-PERMEABILITY 
MATERIAL 

A validated theory for predicting fluxes and forces in a suspension system 
based on closed magnetic circuits employing high-permeability materials was 
reported in Reference 4. The emphasis in Reference 4 was the experimental 
validation of the theory and a subsequent optimization of the bearing surface 
geometry (tooth geometry) to obtain a maximum stabilizing transverse force. 

The purpose here is to present some analytical results that were obtained 
using the above theory for two different tooth configurations: straight 
tooth and circular tooth. The results for these two cases will be presented 
in two separate sections after a brief description of the magnetic quantities 
used. The results include plots of nondimensional values of flux and normal 
and transverse forces as functions of either gap or displacement, with tooth 
geometry (gap * tooth spacing) as the parameter. 

2.1 Description of Magnetic Quantities 

A schematic of the geometry of the teeth on the two bearing surfaces used in 
the analysis is shown in Figure 2-1. It has been shown previously [4] that, 
for most practical situations, the semi-infinite tooth-width approximation 
shown in Figure 2-2 is adequate for the calculation of the fringing magnetic 
field around the edge of the teeth and the normal and transverse forces. 

The limitation of the above analysis, however, is that the transverse displace¬ 
ment in relation to the tooth width and the air-gap thickness cannot exceed a 
certain value > 5) for validity. This, however, was found not to re¬ 

strict too severely the applicability of the theory for practical bearing 
calculations. 

The actual theory and the necessary analytical treatment for the computation 
of fluxes and forces has been described in detail in Reference 4. Only the 
different parameters and variables used in the presentation of the results 
are described here. 
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The two geometrical parameters that appear with the results are and —. It 
is required that the parameter ^ be greater than 0.5 for accurate theoretical 
prediction. The results are reported at various values of the parameter —. 
The main variable is the nondimensional value of the transverse displacement 



The three main quantities computed analytically are <*>, F and F t - These 
given by: 
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( 2 - 1 ) 


where p ± (^) is the permeance of the air-gap per pair of teeth and <t> i 
(Wb) is the flux per tooth. A F (AT) is the mmf present between the 
two surfaces. 


(2) F 


2h F if 
U o (AF) 2 £ 


( 2 - 2 ) 


where F if (N) is the normal force per tooth. 


(2-3) 


The prefix A is used to denote changes in quantities with respect to 
their values in the aligned position. The quantity i appearing in 
Equations (2-1), (2-2), and (2-3) represent the length of the teeth 
for the straight-tooth configuration. 


2.2 Results for Straight-Tooth Configuration and Discussion 

Figure 2-3 represents the variation of air-gap permeance with gap thickness 
for the aligned position of the teeth. If the flux for a given case is known, 
the above curve can be used to determine the corresponding mmf acting across 
the air gap. The variation of normal force with gap is represented in Figure 
2-4 for the aligned-tooth configuration. It is to be noted that the shape of 
tne curve in the above figure is as shown because the ordinate contains the 


parameter 
















With Gap Tliioki 





Figure 2-5 shows the decrease in air-gap permeance from the aligned position 
when the teeth are transversely displaced. The previous figures contain curves 
for different values of the parameter -j|. The drop in air-gap permeance can 
be significant at large displacements and small gaps. 

The decrease in normal force corresponding to the decrease in air-gap per¬ 
meance as the teeth are transversely displaced is shown in Figure 2-6, which 
contains curves for different values of the parameter The normal force 
reduces significantly at large displacement. 

Figure 2-7 shows a plot of the nondimensional value of the transverse force 
with respect to the displacement ratio Curves for different values of the 
parameter — are given in the above figure. The slope of the curves at any 
displacement represents the transverse stiffness. The transverse stiffness 
can be negative at large displacements; this is the result of the interaction 
of flux for the neighboring teeth. 

2.3 Results for Circular-Tooth Configuration and Liscussion 

The analysis procedure for the calculation of fluxes and forces for the 
circular-tooth configuration is the same as that for the straight tooth. In 
all the nondimensional groups mentioned previously the length of the teeth £ 
is replaced by the mean circumference nD where D is the mean diameter. In 
addition, the suffix R is used to denote the circular-tooth case. 

When circular teeth are transversely displaced by an amount e, the actual 
displacement between corresponding segmental elements of the teeth varies 
sinusoidally according to 

e' = e sin6, 

where 9 refers to the angular position. 

In order to compute the fluxes and forces for the circular teeth, the results 
for misaligned, segmental elements (approximated by straight configuration) 
are integrated around the mean circumference. The results presented in this 
section are the integrated values. For the circular tooth, the decrease in 
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air-gap permeance, Che decrease in normal force, and the transverse force v< 
ation with respect to transverse displacement are shown in Figures 2-8, 2-‘ 
and 2-10, respectively. These correspond to Figures 2-5, 2-6, and 2-7 for 
the straight tooth, respectively. Each of the above figures contains curve 
for different values of the parameter When i is replaced by nD the Fig¬ 
ures 2-3 and 2-4 can be used to predict flux and normal force (at e = o) f< 
the circular-tooth configuration. 





































3.0 DYNAMIC MEASUREMENTS FOR CLOSED MAGNETIC CIRCUIT WITH HIG H-PERMEABILITY 
MATERIAL 

A validated theory for predicting fluxes and forces in a closed magnetic 
circuit was reported earlier [4]; the analytical results obtained using this 
theory were presented in Section 2.0 in the form of curves of nondimensional 
values of fluxes and forces versus gap thickness or transverse displacement 
at various values of tooth-geometry parameter. However, the validation was 
performed through measurements under static conditions. 

The purpose of the research described in this section was to detect any pos¬ 
sible effects of transverse vibration and to determine any damping that may 
be present in the system. 

Dynamic measurement of forces was therefore performed on the same magnetic 
circuit used earlier [4]. Comparisons could thus be readily made. Transverse 
oscillations were introduced to the lower core piece (constituting one-half 
the magnetic circuit) using an electronic vibration system. A separate test 
rig was built for this study; this section describes the test rig, the 
associated instrumentation, and the results of the dynamic measurements 
performed. 


3.1 Experimental Apparatus 

3.1.1 Magnetic Circuit 

The magnetic circuit is the same as studied previously [4] and is shown in 
Figure 3-1. A new test rig was constructed to permit the lower core piece to 
be transversely displaced in the Y-direction. A sinusoidal transverse dis¬ 
placement was introduced, using an electrodynamic shaker table. The upper 
core piece,shown in the figure, is fixed to the support housing through two 
three-directional force transducers. The test rig facilitates raising and 
lowering of the upper core piece for varying the gap thickness between the 
two core pieces. 

The magnetizing coil consists of four separate 750-turn coils connected as 
shown in Figure 3-2. This arrangement was implemented in order to allow 
visualization and adjustment of the gap with the coils in place. A total 
current of 1 amp produces 1,500 ampere-turns of magnetomotive force. 






















The core pieces are the same as those used in the previous program. A photo¬ 
graph of one of the core pieces is shown in Figure 3-3. The teeth are 2.5 
inches long, 0.1 inch deep with a center-to-center spacing of 0.060 inch. 

Each pole face contains 28 teeth, and all measurements were performed with 
.020-inch-wide teeth. 

3.1.2 Test Rig 

The test rig is shown schematically in Figure 3-4, and photographs of the rig 
are shown in Figures 3-5 and 3-6. The complete assembly is mounted on a ma¬ 
chine tool bedplate (partially visible in Figure 3-5), which acts as a seismic 
mass. Bearing support for the driven lower core piece is provided by four 
Thomson ball bushing, and pillow blocks which run on two parallel shafts 
attached to the bottom support plate. A schematic of the arrangement is shown 
in Figure 3-7. The pillow-block housings act as clamps which can be adjusted 
to radially preload the bushings. 

The upper core piece is attached to the top mounting plate through two force 
transducers, and this subassembly is supported by four ball-screw/ball-nut 
assemblies which are used to facilitate precise adjustment of the gap height 
between pole faces. Self-aligning pillow-block bearings attach the ball screws 
to the top and bottom support plates. The upper core piece can be positioned 
by simultaneously rotating the four ball screws through the timing belt- 
pulley arrangement shown. The ball screws and the ball bushings are rated 
for 1,000 lbf in the Z-direction and, therefore, were adequate for present 
testing. 

Extreme care was taken to obtain precise parallelism of the magnetic faces 
under both static and dynamic conditions. The top of the bottom mounting 
plate was used as a reference surface for leveling the lower core piece. 
Parallelism between the ball-bushing axes and the reference surface was 
achieved by using shims at the pillow block and shaft mounts. Parallel 
motion of the reference surface along the ball-bushing shafts was verified 
with dial indicators. Following this, the upper core piece and ball-screw 
assemblies were installed with precision gage blocks separating the pole 
faces throughout the assembly operation to assure parallelism. The key step 
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in the assembly for assuring pole-face parallelism was to locate and drill the 
bolt holes for attaching the four ball nuts to the top mounting plate, with 
the gage blocks separating the pole faces. This nullified the potential 
errors in face parallelism due to tolerance stackup of the various machined 
parts. 

To minimize the effects of any external magnetic fields, the top and bottom 
mounting plates were made of nonmagnetic (300 series) stainless steel, and 
the four support plates were made of aluminum. 

The bottom mounting plate was connected to the shaker table by a hexagonal 
steel rod with opposite hand threads at the two ends. This was done so that 
the lower core piece could be displaced horizontally for static measurements 
by simply rotating the rod. Ling Electronics Vibration System was used to 
oscillate the lower core piece. The system and its test cell are shown in 
Figure 3-8. The system could be programmed to provide any arbitrary, periodic 
wave form - including random vibration. For this application, it was used to 
provide a purely sinusoidal excitation at a specified frequency and amplitude. 

3.1.3 Measuring System 

The instruments used for the static measurements and their calibration pro¬ 
cedures are essentially the same as those used in the previous program, for 
which detailed descriptions can be found in Reference 4. However, this sec¬ 
tion briefly describes the instrumentation used for the static measurements. 
Special wave-form analysis instruments used for the dynamic measurements are 
described in Section 4.2. 

Force was measured using two Kistler Type 9251 piezoelectric force transducers 
with a Kistler Model 503 charge amplifier. This system was used also for the 
dynamic force measurements. 

Magnetic flux was measured with an LDJ Electronics Model 701 voltage integrat¬ 
ing flux meter. The relevant static flux measurements were taken from the 
search coil wrapped around the teeth (10 turns), as shown in Figure 3-1. The 
search coil wrapped around the body of the upper core piece was used for 





3-10 




















residual flux measurements while degaussing the core. The procedure used for 
degaussing the core is the same as that described in Reference 4. 

The current source used to energize the magnetizing coils was a Lambda Model 
LL902-OV regulated power supply. During all measurements, the coil current 
was monitored on a digital ammeter, and the power supply regulation was such 
that current could be held constant to within 0.005 amp. When X-Y plots of 
force and flux versus current were generated, the plotter input representing 
current consisted of the voltage drop across a 0.1-ohm precision resistor 
connected in series with the power supply. 

A Schaevitz Model 100 DCD LVDT (Linear Variable Differential Transformer) was 
used to measure the horizontal (tangential) displacement of the lower core 
piece under both static and dynamic conditions. The LVDT is linear over a 
wider range than the capacitance probe system used previously. Figure 3-7 
shows the setup for the LVDT. Its output sensitivity is 101.562 volts per 
inch with a range of ±0.1 inches. Linearity is 0.064 percent of full range. 
The LVDT operates on a ±15-volt DC power supply. 

For all measurements, the pole-face gap was set initially using shim stock 
and then monitored during the test using a .050-inch-range capacitance probe 
with a Wayne-Kerr Model DM100 distance meter. Resolution of this system was 
better than 0.5 percent of full scale, or 250 microinches. The probe was 
attached to an aluminum bracket on the bottom mounting plate, and in this 
manner, any changes in gap due to deflections of the ball screws or ball 
bushings under high magnetic forces were detectable. 

The LVDT and capacitance probe outputs were monitored during all measurements 
using two Data Precision Model 2440 multimeters. 

The data were plotted directly using a Houston Omnigraphic Series 2000 X-Y 
plotter. Dynamic wave forms were monitored on Tektronix Type 504 dual-trace 
oscilloscopes and photographed at relevant data points. 







3.2 Static Measurements 

In the previous research program [4], excellent agreement was obtained between 
experimental measurements and theoretical predictions of magnetic forces and 
flux. These measurements were highly repeatable, giving a high level of con¬ 
fidence in the instrumentation and experimental procedures utilized. 

Following the assembly of the new test rig for the present program, these 
static measurements were carefully repeated as a check on the integrity of 
the new rig. Results were highly affirmative, with the new test rig either 
reproducing the previously measured data or, in some cases, giving improved 
agreement with theory over the previous data. 

However, two sources of error which had to be accounted for were present in 
the measurement of transverse forces. The radial compliance in the ball¬ 
bushing pillow blocks and axial compliance in the ball nuts caused the gap to 
decrease as the current, and thus axial force was increased. At the same 
time, radial compliance in the ball nuts permitted the onset of tangential 
force to slightly shift the horizontal position of the upper core piece. 

Only under conditions when the normal forces were very high (« 300 lbf), the 
data were not perfectly reproducible. At other conditions, the horizontal 
deflection of the upper core assembly was negligible, and the decrease in gap 
could be compensated for by making the gap setting equal to the desired gap 
plus the expected vertical deflections of the ball nuts and ball bushings. 

3.3 Flux and Transverse Force Measurements Under Dynamic Conditions 

One of the primary goals in performing dynamic measurements of the tangential 
force was to quantitatively determine the degree of damping present in a mag¬ 
netic suspension under a given set of conditions. It is hypothesized that the 
tangential or passive transmission of force across the pole-face gap takes 
place through an effective spring and dashpot acting in parallel so that the 
restorative and dissipative force components are directly additive. The 
restorative or spring force is due to the tendency of the field to force the 
teeth into alignment, and the dissipative or damping force is due to eddy 
current and magnetic hysterisis losses which arise from a time varying flux. 
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For an arbitrary, nonlinear force wave form measured at the upper core piece 
load cells in response to a purely sinusoidal oscillation of the lower core 
piece, the energy dissipated per cycle is given exactly by the expression 
(because of the orthoganality relations for the Fourier eigenfunctions): 

E D = n5F o sinQl * (3-1) 

where 6 is the amplitude of oscillation of the lower core piece, F q is the 
amplitude of the first harmonic of the force wave form, and a is the phase 
angle by which the first harmonic of the force wave form leads the purely 
sinusoidal displacement wave form, a is a good indicator of the degree of 
damping in the system. In fact, for a linear system 



where D is the damping coefficient in the system and D £ is the critical damp¬ 
ing coefficient. 

3.3.1 Instrumentation for Wave-Form Analysis 

Measurement of the amplitude of the first harmonic of the tangential force 
wave form F q and its phase angle, with respect to the lower core piece motion 
wave form a, was accomplished using VIC (Vibration Instruments Co.) Model 
235DS Dual-Channel Tracking Analyzer, and a Dranetz Model 305 Phasemeter. 

The tracking analyzer is an electronic system which acts as a narrow band¬ 
width filter, of which the center frequency is precisely and automatically 
tuned to follow (track) the fundamental frequency of a varying, external 
reference signal. The dual-channel unit used here consists of two parallel 
filters tuned to a common reference. The two filter channels are phase 
matched to within 1 degree so the output signals, which are purely sinusoidal, 
can be fed into a phasemeter to measure the phase angle between the input 
signals. 

For the tangential force measurement, the lower core piece motion wave form 
is fed to the reference input and one of the two filter channels (Channel A). 
The tangential force wave form is fed to the second filter channel (Channel B). 
This arrangement is shown schematically in Figure 3-9. The voltmeter provides 
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a digital display of the rtns value of the first harmonic of the force wave 
form F , and the phasemeter gives a digital display of the phase angle a. 


3.3.2 Dynamic Rig Checkout 

3.3.2.1 Tracking Analyzer and Phasemeter . Dynamic data were desired over a 
frequency range of 20 Hz to 100 Hz, so the instrumentation was checked out 
over a much wider frequency range. In setting up the tracking analyzer and 
phasemeter, extensive care was taken to ensure that there would be no super¬ 
fluous indication of phase angle. 

It was determined by an independent check-out procedure that both analyzer 
channels must have the same range setting and that the amplitudes of the 
inputs to the tracking analyzer have the same magnitude. (A resistive 
attenuator can be used to accomplish this.) Under these conditions, the 
phasemeter could be zeroed to ±0.2 degrees. 

3.3.2.2 Complete Rig Checkout . Check-out tests were performed with a 10-mil 
gap at flux levels of 90 and 135 kilolines. Lower core piece oscillation 
frequencies ranged from 20 to 80 Hz, with amplitudes ranging from ±.004 to 
±.010 inch. 

Initially, structural vibration problems were encountered at frequencies above 
50 Hz, due primarily to radial (horizontal) compliance in the ball-screw/ball- 
nut assemblies. The excitation was through the magnetic force transmission 
at the pole-face gap and a small dynamic force transmitted radially through 
the ball bushings due to the slight error in colinearity between the ball¬ 
bushing axes and the shaker-drive axis. This structural, vibration-threshold 
frequency was increased from 50 Hz to 70-80 Hz, giving a range of 20-60 Hz 
over which dynamic data could be acquired with a high level of confidence by 
implementing clamping bolts between the side supports and the top mounting 
plate. 


In the initial attempt at measuring the phase angle between tangential force 
and displacement, the LVDT was used to generate the displacement signal, and 
the indicated phase angles appeared suspiciously large. However, by measuring 







the phase angle between the LVDT signal and the signal given by an accelero¬ 
meter mounted on the lower core piece, it was established that the major 
portion of the indicated phase angle between force and displacement was due 
to differences in the electronic phase shifts introduced by the LVDT and the 
charge amplifier used with both the accelerometer and the force transducers. 

This situation was remedied by basing the phase measurement on a comparison 
between the force transducer and accelerometer signals by using the LVDT as 
an indication of motion amplitude only. Kistler Model No. >03 charge ampli¬ 
fiers were used with both the force transducers and the accelerometer. The 
two charge amps were found to produce zero net phase shift as long as their 
range and sensitivity settings were identical. This was established by pair¬ 
ing the two identical force transducers with two charge amps and measuring 
the phase angle between the charge-amp outputs. While the raw transducer 
signals are naturally in-phase since both transducers are mounted between the 
same rigid members, it was possible to generate a significant indicated phase 
angle between the charge-amp outputs by varying the sensitivity of one charge 
amp relative to the other. With identical sensitivites, the phasemeter indi¬ 
cated zero to within ±1° at frequencies up to 70 Hz, giving the complete 
phase-measuring system a resolution of ±1°. 

3.3.3 Experimental Procedure 

The connecting rod between the shaker and the lower core piece support was 
adjusted so that the oscillation took place about the position of complete 
tooth alignment. The static LVDT output level corresponding to this position 
was marked on an oscilloscope, and with the rig operating, the connecting rod 
was adjusted for each data record until the DC level of the LVDT-output wave 
form coincided with the marking. 

The midpoint of the shaker-table motion showed a tendency to drift somewhat, 
particularly at lower frequencies and amplitudes. This drift tended to be 
unpredictable and could not always be compensated for by connecting-rod ad¬ 
justments, resulting in some inaccuracies in the data. Furthermore, while 
the frequency of the shaker-table motion showed no tendency whatever to 
drift from the set value, the motion amplitude did tend to drift somewhat, 
making it necessary to record the LVDT-output amplitude. 

TCCHMOtOO* 
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An additional source of error was due to the variation in axial force with 
tangential position and the resulting gap variation. 

3.3.4 Experimental Results and Discussion 

3.3.4.1 Phase Angle . Phase angle between the tangential force and the tan¬ 
gential displacement was measured under the following conditions: .005- to 
.020-inch gap, 0 to 60-Hz frequency, 0 to + .010-inch amplitude of transverse 
oscillation, and at 90- and 135-kilolines flux. In all of the above cases, 
the phase angle measured was less than or equal to the resolution of the 
phase angle measuring system (resolution is approximately +1°). Under these 
conditions, even if the system possessed some damping (which is within the 
limit of resolution of the measuring system), the magnitude is extremely 
small and inadequate for any practical bearing application. 

3.3.4.2 Tangential Force . Plots of the measured, tangential-force amplitudes 
as a function of the measured motion amplitudes are shown in Figures 3-10 and 
3-11. The corresponding measured static forces are shown on each plot for 
comparison. It is notable that in most of the cases shown there is a marked 
decrease in the force-amplitude level with increasing frequency at a fixed- 
motion amplitude with the static case corresponding to zero frequency. Even 
though there are sources of inaccuracy in the previously described data which 
are difficult to quantify, the data appear to show with sufficient regularity, 
particularly the cases of 90 kilolines at .005-inch gap and 135 kilolines at 

-inch gap, that the phenomenon of a decrease in the passive force response 
with increasing frequency is real. 

Figure 3-12 shows the tangential force wave form at a motion amplitude of 
.005 inch (a normal motion level) and an amplitude of .016 inch (an extreme 
level in view of the .020-inch tooth width). The apparent, pure sinusoidal 
character of the .016-inch wave form confirms the predominant linearity of 
the tangential force-displacement characteristic over a wide range of motion 
amplitudes. 
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4.0 ANALYSIS OF PERMANENT MAGNETS IN REPULSION 




The previous sections of the report described the study performed on closed- 
circuit, electrically energized magnetic systems employing high-permeability 
material. 

This and the following sections of the report describe the study of permanent 
magnets in both open- and closed-circuit systems. 

With the progress in the science and technology of magnetic materials, perma¬ 
nent magnets are being utilized in a variety of applications, some of which 
involve heavy loads [5]. If constant field strength of high values is the 
virtue of the permanent magnet, electromagnet has the advantage of ease 
with which its field strength can be varied (simply by varying the current 
through the coil). In any practical application, the two kinds of mangets 
can be used in combination to derive the advantages of both. 

In this section, a report is made of the analytical study performed on open- 
circuit magnetic suspensions. The objective here is to calculate both the 
normal and transverse forces between two repelling permanent magnets. Because 
of their significantly enhanced properties, high-coercivity rare-earth-cobalt 
permanent magnets are considered here. Although the approach is fairly gen¬ 
eral, the analysis procedure discussed is more convenient for the analysis of 
high-coercivity permanent magnets. 

It has been shown by Tsui et al [6] that, for the analysis of high-coercivity 
permanent magnets (like those of samarium cobalt or praseodymium cobalt) in 
repulsion, the intrinsic material properties are the relevant ones. In re¬ 
pulsion applications the magnetization or the intrinsic induction M is the 
important property. For the rare-earth-cobalt alloys, Senno and Tawara [7] 
have shown that the intrinsic hysterisis loops (M versus H) are very nearly 
square in shape. This means that the magnetization M is a constant for values 
of the demagnetizing field, less than the coercive field intensity. The 
relation between the flux density B and the field intensity H for a rare- 
earth-cobalt magnet, therefore, takes on the form 
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B = H H + M, 


(4-1) 
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where M is independent of B and H, and p is the permeability of the surround¬ 
ing space. The relative permeability for the magnet has been assumed to be 
unity. 

In magnetostatic field problems it is customary to use a vector potential 
from which the field parameters can be calculated. Thus the vector potential 
A may be defined as [1] 


V x A = B 


(4-2) 


where V is the vector differential operator. The requirement that the diver¬ 
gence of vector B be zero (one of Maxwell's electromagnetic field equations) 
is automatically satisfied by Equation (4-2) because of a vector identity. 

In a region of uniform permeability if I is the current density, it has been 
shown that [1] 

<«-» 

where p v is the permeability and r is the magnitude of the radius vecotr r. 

In permanent magnets, the magnetization M can be thought of as the magnetic 
moment per unit volume of the permanent circulating currents. The vector 
magnetic potential for M is therefore written as [1]: 
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(4-4) 


By using vector identities and Gauss's Theorem, Equation (4-4) can be written 
as [1] 


i H 5 *** / H*. 


(4-5) 


where n is the unit vector normal to the surface S bounding the volume V. 

The first integral on the right-hand side of Equation (4-5) is a volume in¬ 
tegral and therefore, by comparing this with the integral in Equation (4-3), 
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V x M can be thought of as the volume current density. V x M has to be 

2 

divied by in order to get the familar units of (A/m ). Therefore, the 
equivalent volume-current density for a permanent magnet due to M is 


J 


(4-6) 


_ 2 

where J is in Amp/(meter) . Advancing a similar treatment, an equivalent 
surface-current density for a permanent magnet due to M is 


K 


M x n 
y v ' 


(4-7) 


where K is in Amp/meter. 


The magnetization M in a permanent magnet can therefore be represented equiva¬ 
lently by a volume-current density given by Equation (4-6) and a surface- 
current density or current-sheet density given by Equation (4-7). This kind 
of representation for M is amenable to mathematical analysis. 

It has been assumed for purposes of analysis performed here that the magnet¬ 
ization M is uniform (i.e., unidirectional) and of constant magnitude through¬ 
out the volume of the magnet. This assumption appears very reasonable for 
the rare-earth-cobalt magnets. The nonuniformity is partly introduced during 
the process of compacting the powder alloy, and therefore the extent of non¬ 
uniformity depends on the size of the magnets. However, this nonuniformity 
of magnetization throughout the volume of the magnet was quoted by the manu¬ 
facturers to be less than 5 percent. With the assumption of uniform magnet¬ 
ization, it can be easily seen that the volume-current density becomes zero 
and the surface-current density assumes a constant magnitude over the entire 
magnet surface. 

4.1 Analysis Procedure 

Under the assumption of uniform magnetization, the equivalent surface currents 
enveloping the nonpolar surfaces of the permanent magnet form a current sheet 
of constant density. This current sheet can be considered as the summation 
of a series of electric circuits. Such electric circuits are present on the 
nonpolar surfaces of every magnet in the system. 
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With such a representation of the permanent magnet, calculation of forces 
basically involves the use of Ampere's Force Law [1], which predicts the 
force between two electric circuits. The force is given by the integral of 
the double cross product: 


F = u o ll’ <j) $ 


ds x (ds 1 x r) 
4ir r3 


(4-8) 


where ds and ds' are the elemental arcs of the two electric circuits and r is 
the radius vector between ds and ds'. I and I' are the currents flowing in 
the two circuits. Using a vector identity, Equation (4-8) can be written as: 


F = - 


Uq 11 ' 
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r (ds . ds') 
r3 


(4-9) 


u , which is the permeability of the surrounding space, is used in Equations 
(4-8) and (4-9) because the relative permeability for rare-earth-cobalt magnets, 
which is very nearly unity, has been assumed to be 1.0. It can be noticed 
from Equation (4-9) that, when ds and ds' are perpendicular to each other, no 
force exists between the two circuit elements. 


The analysis procedure, therefore, involves a repeated use of Equation (4-8). 
The sign of each integral has to be properly considered. This approach to 
repulsion force calculation has been used earlier by Borcherts [8] and others 
16, 91. 


4.2 Analysis for Rectangular Parallelopiped-Shaped Magnets 

Figure 4-1 shows two rectangular parallelopiped-shaped permanent magnets in 
a repulsive mode. The two magnets are assumed identical in shape in order to 
reduce the computational effort. The linear dimensions are nondimensionalized 
with respect to the width of the magnet. The two magnets are shown misaligned 
by an amount E in the X-direction. Study of misalignment is of importance 
in bearing applications because of the associated changes in forces and 
stiffness coefficients which determine the stability of the suspension system. 
In order that they produce repulsion, the directions of magnetization M of the 
two magnets are taken to be in +2 and -Z directions. 















The equivalent surface currents due to magnetization M' of the upper magnet 
are shown as K' per unit surface area of the magnet. The direction of these 
currents is as shown in Figure 4-1. Similarly, the equivalent surface cur¬ 
rents K per unit surface area of the lower magnet are as shown in the figure. 
The points on the upper magnet are denoted by (X', Y', Z’), while (X, Y, Z) 
represents points on the lower magnet. 


If we now consider elements on the surfaces of each of the two magnets, the 
repulsive or attractive forces between any two elements can be calculated by 
using Ampere's Inverse Square Law, Equation (4-9). 


Summing all the elemental forces over the faces 1 and 1', the repulsive force 
Fj_^, is obtained as: 
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(4-10) 


where r^_^, is the position vector between the points (X, Y, Z) and (X', Y', 
Z'). Therefore, r^_ lf is given by: 

*1-1' = < X, - X >*x + 0 ’ ° y + <*'-Z>S B . < 4 - U > 

where U , U , and U are the unit vectors in X-, Y-, and Z-directions, respect- 
x’ y’ z _ 2 

ively. The magnitude of the vector Is rj_^,. In Equation (4-10), W 

appears outside the integral sign because of nondimensionalization. 

By symmetry, the repulsion force between faces 3 and 3' F^^i is equal to F^_^,. 
Similarly, 
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and 


E U x + (Y’-Y) U y + (Z'-Z) U z . 


(4-13) 


Again by symmetry, F 2 _ 2 , = F 4 _ 4 ,. The attraction force F^, is given by: 
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J l-3' = (X '“ X) C x + ™ y + (Z '“ Z) U z ; 
F l- 3 ' is equa l to F^i _3 by symmetry. 

The attraction force F 2 _^, is given by 
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The attraction force 
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(4-19) 


A vector sum of Z-components of the above forces yields the normal component 
of the net repulsive force. The transverse component of the net repulsive 
force is obtained by summing vectorially the X-components of the above forces. 
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All the above quadruple integrals can be integrated analytically twice to re¬ 
duce them to double integrals in Z and Z'. The normal component of the total 
repulsive force F z can be written as 

U 2 0 L+EG 

F z = 7- KK’W / / f (E, T, Z'-Z) dZ’ dZ (4-20) 

-L EG 

The variables Z and Z' in Equation (4-20) can be transformed to 4 and n by 
4 = Z'-Z 
n = Z'+Z. 

Since the integrand in Equation (4-20) is only a function of (Z'-Z), integra¬ 
tion with respect to n after transformation becomes trivial and, therefore, 
the double integral of Equation (4-20) can be reduced to a single integral 
given by 
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Similarly, the transverse component of the total repulsive force F is given 
by 
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The functions f z and f^ are given by the following: 
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The above integrals for the normal 


numerically. The details of the numerical 


nd transverse forces were evaluated 


evaluation will be given in a later 
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4.3 Analysis for Cylinder-Shaped Magnets 

Figure 4-2 shows two cylinder-shaped, permanent magnets in a repulsive mode. 
The two magnets are assumed to have identical dimensions. All the linear 
dimensions are nondimensionalized with respect to the radius R. The magnet¬ 
ization vectors M and M' are in the +Z and -Z directions as shown in the 
figure. Equivalent surface currents K and K' are also shown in the figure. 
As in the previous section, (X, Y, Z) and (X', Y', Z') represent points on 
the lower and upper magnets, respectively. 


Considering elements on the surfaces of each of the magnets, the repulsive 
force between any two elements can be calculated by using Ampere's Law, 
Equation (4-9). A summation of these elemental forces gives the total repul¬ 
sive force: 

p ? 0 L+EG 2tt 2 it 

F = KK'R / f f f ~3 cos((t>-<*>') d<(> 'd<}> dz'dz (4-23) 
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and 

r = (E + cos<f>'-cos<t>)U x + (sin(|> '-sin<J>)U y + (Z'-Z)U z . (4-24) 


The X-component of F gives the transverse force while the Z-component yields 
the normal force. By symmetry, the Y-component is zero. The above integral 
can be integrated analytically twice with respect to Z and Z', yielding: 
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During the process of integration, when E=0, the parameter P can take on a 
value of zero. Under these conditions, f z and f^ as defined above will be 
infinite. In order to overcome this problem, the above integrals were 
evaluated separately for P=0 and are given by: 
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where 


P* = E 2 - 4E sin2X sin2Y + 4sin 2 2Y 

P, Q, and are as defined earlier, but X and Y substituted in place of <f> and 
if’ using Equation (4-29). 

Similarly, 

u tt/ 2 ^ n-Y 

F x = 4^ - 8 / dY / lf x (P) + f x (P*)J dX 

0 Y 

where 

P* = E 2 - 4E sin2X sin2Y + 4 sin 2 2Y. 

P, A, and f^ are as defined earlier. 


Computational details for the above integrals will be given in a later section. 


4.4 Results and Discussion for Rectangular Parallelopiped-Shaped Magnets 

The integrals for the normal and transverse forces given by Equations (4-21) 
and (4-22) were evaluated numerically by using a 10-point Gaussian quadrature. 

A computer program was written in order to perform this numerical evaluation. 

The values computed were the nondimensional values of the normal and transverse 
F z F x 

forces - and -, respectively. 
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The computer program was run on the General Electric Mark IV computing 
facilities available at MTI. The numerical technique used (i.e., the 10- 
point Gaussian quadrature) was accurate to within 0.01 percent. 

Using the above computer program, a systematic, parametric study was performed. 
The parameters studied were the eccentricity E, the gap EG, the thickness T, 
and tiie length L. All these parameters have been nondimensionalized with re¬ 
spect to the width (W) of the magnet. 










The calculated normal and transverse forces are plotted against the operating 
parameters E and EG. T=1.0 and 5.0 and L=.2, .5 and 1.0 are studied which 
represent a practical range of these geometric parameters. Figures 4-3 through 
4-26 show the results for the case of two rectangular parallelopiped-shaped 
magnets. 

Figure 4-3 shows the normal component of the repulsive force plotted against 
the eccentricity E. The positive value of Z-force is repulsion and the nega¬ 
tive value is attraction. For any constant value of gap EG, the force is seen 
to decrease with eccentricity, as expected. Also, for values of E greater 
than about 0.85, the force changes sign. This is because, at large eccentri¬ 
cities, faces of opposite polarity come closer to each other than the distance 
between faces of like polarity. This attraction force at large eccentricities, 
however, becomes insignificant at large values of gap. 

Figures 4-4 and 4-5 show Z-force (the normal component of the repulsive force) 
plotted against the eccentricity E. Figures 4-3, 4-4, and 4-5 have the same 
value of T, but increasing values of L; therefore, the cross-sectional area 
of the magnets, remains the same. It can be seen from these figures that 
the Z-force (for example at EG=0) does not increase proportionately with in¬ 
crease in length. If we consider that a longer magnet is made up from stack¬ 
ing shorter magnets one on top of the other, the mutual demagnetizing effects 
are evident. Therefore, the Z-force increases with increase in L, but less 
than proportionately. 

Figures 4-6, 4-7, and 4-8 show plots of Z-force versus eccentricity at T=5.0 
and various values of L. If we compare figures 4-6, 4-7, and 4-8 with Figures 
4-3, 4-4, and 4-5, the effect of increasing the thickness can be observed. 
Increasing the thickness while keeping the length the same is equivalent to 
stacking the magnets one by the side of another. The mutual interaction and 
demagnetizing effects are clearly evident. For this reason, the Z-force in¬ 
creases with T, but less than proportionately. By comparing the results shown 
in Figures 4-3 through 4-8, it can be seen that, for the same volume of the 
magnet, increasing the thickness increases Z-force more than increasing the 
length. This is because the effective gap for a longer magnet is larger than 
that for a magnet with higher value of T. 




























































The transverse component of the repulsive force versus eccentricity E at 
various values of the gap is shown in Figure 4-9. The transverse force is 
zero at E = 0 because of symmetry. But, interestingly enough, the transverse 
force reaches a maximum and then starts decreasing with further increase in 
eccentricity. The transverse component of the repulsive force depends on two 
quantities: the angle (0) between the total repulsive force and the X-axis, 
and the distance (r). r increases with eccentricity while 0 decreases. The 
repulsive force is inversely proportional to the square of the distance r, 
whereas the transverse component is proportional to the cosine of the angle 0 
Because of these two competing quantities, the X-force reaches a maximum at a 
certain value of E and then starts decreasing. The value of E for which the 
X-force is a maximum depends on the gap, and this value of E increases with th 
gap, as can be expected. 

At small values of eccentricity E, the slope of the repulsive X-force versus 
curve E is positive. In this region, the equilibrium is unstable. However, 
it can be seen from Figures 4-9, 4-10, and 4-11 that, at far enough eccentri¬ 
cities, the equilibrium can indeed be stable. Although a stable equilibrium 
can be established in the X-and Z-directions, instability cannot be avoided 
in the Y-direction (perpendicular to X-Z plane), simultaneously. 


Figures 4-10 and 4-11 show plots of X-force versus eccentricity at various 
values of gap for T=1.0 and L=.5 and 1.0, respectively. Figures 4-12, 4-13, 
and 4-14 show similar plots but for T=5.0. An observation of these figures 
reveals the following: (1) at a given value of T, an increase in length pro¬ 
duces transverse forces which increase more than proportionately to begin 
with; (2) increase in the thickness T produces a more than proportionate in¬ 
crease in the transverse force. These observations can be rationalized by 
considering that the larger magnets are made by stacking smaller magnets. 

It can then be easily seen that the mutual interaction of smaller magnets 
enhances the transverse forces under some conditions. 

Figures 4-1!) through 4-20 show plots of Z-force versus the gap EG at various 
values of the eccentricity E and at different values of L and T. The positive 
values represent repulsion, while negative values represent attraction. The 
Z-force goes to zero at large values of gap, as expected. 
















Fig. 4-10 X-Force Versus E (T = 1.0, L = 0.5) 
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Figures 4-21 through 4-26 show the variation of X-force with the gap EG at 
various values of the eccentricity E and at different values of T and L. The 
transverse forces approach zero at large values of gap, as expected. The rate 
of decrease appears to depend on the size of the magnet. 

4.5 Results and Discussion for Cylinder-Shaped Magnets 

The double integrals for the normal and transverse forces were evaluated by 
using the Gaussian quadrature technique. Considerable computational diffi¬ 
culty was experienced in terms of convergence of the integral at lower order 
Gaussian quadrature schemes. This difficulty was overcome by dividing the 
interval of the outer integral (of the double integral) into four equal sub¬ 
intervals. A 24 x 16 quadrature was used in each of these subintervals. A 
separate computer program written for this purpose was run on a Control Data 
computer system. 

A parametric study similar to the one reported in the previous section was 
carried out for the cylinder-shaped magnets also. The parameters studied 
were E, EG, and L, nondimensionalized with respect to the radius R of the 
magnets. The results have been plotted with respect to E and EG. Figures 
4-27 through 4-38 show these plots. 

Figure 4-27 shows a plot of the normal repulsive force against eccentricity. 

No negative values of Z-force are predicted in this case, as opposed to the 
rectangular magnets. This appears to be due to the geometry of the cylin¬ 
drical magnets. 

Figures 4-28 and 4-29 show the variation of Z-force with eccentricity at 
various values of gap and at different values of L. An Increase in L increases 
the Z-force, but less than proportionately. The reasons are similar to those 
cited in the previous section on rectangular magnets. 

Figure 4-30 shows the transverse component of the repulsive force (X-force) 
plotted against eccentricity at various values of the gap. The transverse 
force is zero at zero eccentricity because of symmetry. The transverse force 
increases with E, reaches a maximum, and then decreases with further increase 
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Fig. 4-27 Z-Force Versus E (L = 0.2 
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in E. Reasons for such a behavior are identical to the ones cited in the pre¬ 
vious section. In the decreasing side of the curve while stability of equili¬ 
brium can be achieved in X- and Z-directions, instability cannot be avoided 
in the Y-direction simultaneously. 

Variations of X-force with eccentricity at values of L = 0.5 and 1.0 are shown 
in Figures 4-31 and 4-32. An increase in L initially produces a more than 
proportionate increase in X-force. By considering that the larger magnets 
are made by stacking smaller magnets one on top of the other, it can be seen 
that mutual interaction may enhance the transverse force under certain 
conditions. 

X-force at various values of L is shown plotted against the gap EG in Figures 
4-33 through 4-35. The repulsive Z-force goes to zero at large values of gap 
for all eccentricities, as expected. A similar behavior of X-force with gap 
is shown in Figures 4-36 through 4-38. 







































































5.0 EXPERIMENTS WITH PERMANENT MAGNETS IN REPULSION 


The analytical calculation of forces between two repelling permanent magnets 
was reported in the previous section. In this section, a report is made of 
the experimental work performed regarding the measurement forces between 
repelling permanent magnets to validate the theory reported in the pre¬ 
vious section. Both normal and transverse forces were measured using 
quartz piezoelectric transducers as functions of gap and transverse dis¬ 
placement. A direct comparison is made of the forces predicted by the 
theory and those measured experimentally. 

5.1 Experimental Apparatus 

5.1.1 Magnet Material 

The permanent magnets that were used in the experimental study were rectangular 
parallelopiped-shaped and were made of high coercivity, rare-earth-cobalt 
(Hicorex 90A from Hitachi Magnetics Corporation). A list of some properties 
of this material appears in Table 5-1. 

5.1.2 Test Rig 

The test rig used is the same as that described in Section 2.0 and is shown 
in Figure 5-1. Each of the two magnets being tested is held in a separate 
aluminum holder as shown in Figure 5-2. Extreme care was taken not to tighten 
the bolts much, because of the extreme brittleness of the sintered, rare-earth- 
cobalt material. The aluminum holders were directly attached to the top 
(movable) and the bottom plates, as can be seen in Figure 5-3. The top plate 
was moved up and down over the screw rods and the ball bushings, thereby 
changing the air gap separating the two magnets. The magnets were transversely 
displaced by manually moving the translation stage, mounted on the bottom 
plate. It was found necessary to mount a translation stage on the bottom 
plate in order to achieve a transverse displacement of 2 inches. 

The air gap and the transverse displacement were monitored by dial indicators 
and also by DC-LVDT's as seen in Figure 5-1. Any required gap was set by us¬ 
ing nonmagnetic shims and the dial indicator. The outputs of the LVDT's were 
used for obtaining plots of forces versus gap or transverse displacement. The 




TABLE 5-1 

PROPERTIES OF RARE-EARTH-COBALT MATERIAL 
(HICOREX 90A OF HITACHI MAGNETICS CORPORATION) 

Typical Room-Temperature Properties 


Residual Induction (B^), gauss 8,200 

Coercive Force (H c >, oersteds 7,500 

Intrinsic Coercive Force oersteds >30,000 

Energy Product (B.H) Max., mgo 16 

Recoil Permeability 1.05 

Required Magnetization Field (H ) oersteds 20,000* 

Hardness, Vickers D.P.H. 500 

Density grams/cm'* 8.2 

Tensile Strength, psi 5,000 

Resistivity, p, micro ohm-cm (@ 25°C) 50 

Thermal Conductivity, K, (cal-cm)/(°C.sec.cm 2 ) 0.025 

Thermal Expansion || to Orientation, cm/cm, °C 5 x 10 ^ 

Thermal Expansion J_ to Orientation, cm/cm, °C 13 x 10 ^ 


*Magnetization requirement based on having a magnet that has been thermally 
demagnetized. Higher fields are required if demagnetized state was achieved 
by a magnetic field. 


Ref.: "Hicorex Rare-Farth-Cobalt Permanent Magnets", Technical 
Brochure-Hitachi Magnetics Corporation, Michigan. 
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magnets were accurately aligned in the transverse direction by searching for 
a minimum value of the transverse force. Figure 5-4 shows a view of the test 
setup with all the instrumentation used for the measurements reported here. 

5.1.3 Force Measurement 

Both the normal and the transverse forces are measured using two multi¬ 
directional quartz piezoelectric force transducers, and the measurement was 
performed in a way similar to that described in Section 3.0. 

The outputs of the two force transducers are connected to the charge amplifiers 
through summing junctions so that the normal and transverse forces are measured 
directly. The outputs of the charge amplifiers are direclty connected to the 
Y-axis of an X-Y plotter, whenever plots are to be recorded. The input cables 
to the charge amplifiers are always kept clean to avoid stray capacitance and 
the resulting drift of the charge amplifier output. 

5.2 Results and Discussion 

Experimental results will be presented in the form of plots of the normal and 
the transverse force versus gap or transverse displacement. 

5.2.1 Z-Force Versus Gap 

Figure 5-5 shows the curve of normal force versus gap when the two magnets 
are perfectly aligned. The figure snows the curves for both decreasing and 
increasing gaps, although the two curves do not significantly deviate from one 
another. 

Figure 5-6 shows curves of Z-force versus gap (similar to that shown in Figure 
5-5) at various values of transverse displacement (e). The minimum gap reached 
in all the tests was .010 inch. The forces reduce significantly at large 
values of transverse displacement and, as can be seen in the figure, the Z- 
force in fact becomes attractive instead of being repulsive at very large 
values of e. 
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5.2.2 X-Force Versus Displacement 

The transverse force versus transverse displacement is plotted in Figure 5-7 
at various values of gap. The forces reduce significantly at large values of 
gap. There appears to be a relatively large range of displacement over which 
the variation in X-force is minimal. The maximum transverse force appears to 
be only about 30 percent of the maximum normal force. 

5.2.3 Z-Force Versus Displacement 

Figure 5-8 shows a plot of Z-force versus transverse displacement at various 
values of gap. Although the curves shown in the figure were independently 
obtained, they are, in effect, cross plots of the data shown in Figure 5-6. 
Figure 5-8, however, clearly shows the variation of normal force as the two 
magnets are displaced transversely, one with respect to the other. At large 
values of displacement, the normal force becomes attractive instead of being 
repulsive. 

5.2.4 Comparison with Theoretical it as 

In order to make a direct comparison between the theoretical predictions and 
the experimental results, the computer code that was written to calculate and 
plot the nondimensional forces versus the nondimensional values of gap and 
transverse displacement (reported in Section 3.0) was run for the values of T 
equal to 1.0625 and L equal to 0.2469, corresponding to the magnets used in 
the experiments. The results thus obtained are shown by solid lines in Fig¬ 
ures 5-9 through 5-11 and are similar to the results shown in Section 3.0. 

In these figures, the vertical axes on the left show the nondimensional values 
of the forces and the vertical axes on the right show the corresponding dimen¬ 
sional values. Both the nondimensional and the actual values of either the 
gap or the displacement are shown on the horizontal axes. 

For comparison purposes, the theory is made to fit the experimental data for 
normal force at E = 0, and EG = .05, and the equivalent surface-current densi¬ 
ties of the magnets (assumed equal to each other) are deduced as indicated in 
Equations (5-1) through (5-4) which follow. 
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Fig. 5-9 Z-Force Versus Cap ( I = 1.0625, L = 0.2469 ) 
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Fig. 5-10 X-Force Versus Displacement (T = 1.0625, L * 0.2469) 






















The nondimensional Z-force is given by 


assuming K = K', 




(5-1) 


(5-2) 


From Figure 5-11, I z = 1.5162 at E = 0 and 
F z = 36.1 lbf for the same values E and EG. 
(5-2) and for p = 1.05*4ir*10 2 ~ ^A/m) t * le 
ties of the magnets are given by 


EG = .05. Figure 5-11 also shows 
Using these values in Equation 
equivalent surface-current densi- 


K = K' 


6.2519 * 10 J . 


(5-3) 


The catalog value of the equivalent surface-current densities can be calculated 
by 


K = K’ 


Mr = Br 
P P 


-7 (Wb/m ) 
(A/m) 


= 6.2146 


(5-4) 


The values of K and K', shown in Equations (5-3) and (5-4), differ by less 
than 0.6 percent, and this demonstrates an excellent agreement between theory 
and experiments. 


p and the associated values of K and K.' can also be calculated so that K and K' 
deduced from force measurements (as done in Equation 5-3), equal to those 
calculated from the catalog value of p r (as done in Equation 5-4). The values 
of p, K and K*calculated following this procedure are given by 


P 


1.04 


4 " 


(Wb/m 2 ) 

(A/m) 


(5-5) 
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and 


K 


6.2744 * 10 meter 


(5-6) 


An excellent agreement between theory and experiments can again be seen, because 
the value of p calculated in Equation (5-5) is very close to the catalog value 
of 


1.05 * 4tt * 10' 


(Wb/m 2 ) 

(A/m) 


Using the experimentally deduced values of K and K' given by Equations (5-3) 
and (5-4), the experimental results described in the previous sections are 
plotted (marked by data points) in Figures 5-9 through 5-1], mentioned 
earlier. The vertical axes on the right show the actual values of the forces. 

A very good agreement between the theoretical predictions and the experimental 
results can be noticed in these figures. In addition to predicting the normal 
forces accurately, the theory is able to predict the transverse forces with 
great accuracy. 


Thus, it can be stated that the current-density analogy for the permanent 
magnets used in theoretical predictions is a valid and usable approach to 
predict forces in open-circuit-type magnetic suspensions employing permanent 
magnets alone. 


Based on the findings described above, it appears that, although nonuniformi¬ 
ties in a magnet do exist, their influence on force predictions in open-circuit, 
permanent-magnet systems is only of secondary importance. 
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6.0 PERMANENT MAGNETS IN CLOSED MAGNETIC CIRCUIT 


In this section, a report is made of the study performed on magnetic suspen¬ 
sions based upon closed magnetic circuits, energized by permanent magnets. 
Because of their high coercivity, rare-earth-cobalt permanent magnets were 
used in the experimental study. 

The theory that was developed for prediciting forces in cl< ;ed-eireuit- 
type magnetic suspensions has been reported earlier. This theory has been 
validated under both static and dynamic (vibratory) conditions, but the magnetic 
circuit considered was energized by electrical means only. However, in order 
to reduce the electrical power requirements, it may be beneficial to energize 
the magnetic circuits by permanent magnets alone or in combination with 
electrical means. The objective of the research reported here has been to 
compare the forces measured experimentally with theoretical predictions for 
closed-circuit-type magnetic suspensions energized by permanent magnets. 

Only plane geometry was considered in the experiments, and a slightly modified 
version of the test apparatus described in Sections 3.0 and 5.0 was used in 
the present study. 

6.1 Magnetic Circui t 

The magnetic circuit studied was essentially the same as that described in 
Reference 4 and referred to in Section 3.0 of this report. The coils used 
for energizing the magnetic circuit were removed and instead, permanent mag¬ 
nets were installed in series with the magnetic circuit as shown in Figure 
6-1. The reluctance of the permanent magnets, therefore, plays an important 
role in determining the circuit performance. Provision was made for instal¬ 
ling up to four magnets of approximately 0.5-inch thickness each. The toler¬ 
ance on the thickness of the magnets was 0.0001 inch and the surface finish 
was better than 16 microinches; these features enabled the magnets to be in¬ 
stalled in the magnetic circuit with a reasonably tight fit, ensuring no air 
gap between the magnets and the pole pieces. 

The core piece with machined teeth used in earlier experiments (.020 inch 
thick and .060 inch pitch) was modified for the present study so that the 
legs of the core piece could be used as pole pieces for the present work. 
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6.0 PERMANENT MAGNETS IN CLOSED MAGNETIC CIRCUIT 


In this section, a report is made of the study performed on magnetic suspen¬ 
sions based upon closed magnetic circuits, energized by permanent magnets. 
Because of their high coercivity, rare-earth-cobalt permanent magnets were 
used in the experimental study. 

The theory that was developed for prediciting forces in clo ;ed-circuit- 
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validated under both static and dynamic (vibratory) conditions, but the magnetic 
circuit considered was energized by electrical means only. However, in order 
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the magnetic circuits by permanent magnets alone or in combination with 
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compare the forces measured experimentally with theoretical predictions for 
closed-circuit-type magnetic suspensions energized by permanent magnets. 

Only plane geometry was considered in the experiments, and a slightly modified 
version of the test apparatus described in Sections 3.0 and 5.0 was used in 
the present study. 

6.1 Magnetic Circuit 

The magnetic circuit studied was essentially the same as that described in 
Reference 4 and referred to in Section 3.0 of this report. The coils used 
for energizing the magnetic circuit were removed and instead, permanent mag¬ 
nets were installed in series with the magnetic circuit as shown in Figure 
6-1. The reluctance of the permanent magnets, therefore, plays an important 
role in determining the circuit performance. Provision was made for instal¬ 
ling up to four magnets of approximately 0.5-inch thickness each. The toler¬ 
ance on the thickness of the magnets was 0.0001 inch and the surface finish 
was better than 16 microinches; these features enabled the magnets to be in¬ 
stalled in the magnetic circuit with a reasonably tight fit, ensuring no air 
gap between the magnets and the pole pieces. 

The core piece with machined teeth used in earlier experiments (.020 inch 
thick and .060 inch pitch) was modified for the present study so that the 
legs of the core piece could be used as pole pieces for the present work. 











During the modification, it was necessary to maintain the alignment of the 
teeth on the two legs of the core piece. In order to accomplish this, 
aluminum plates were fixed to the core piece on either side by means of dowel 
pins and screws, before machining away the material of the core piece. 

Figure 6-2 shows a photograph of the modified core piece with a permanent mag¬ 
net installed. The three pieces seen adjacent to the magnet in Figure 6-2 
are high-permeability iron pole pieces. These are used whenever different 
levels of magnetization have to be realized. The aluminum plates and dowel 
pins on the sides of the core piece can be clearly seen in this figure. 

6.2 Experimental Apparatus and Procedure 

6.2.1 Materials 

The permanent magnets used to energize the magnetic circuit in this work were 
the same as those used for repulsion experiments. The properties of this 
rare-earth-cobalt material have already been listed in Table 5-1 of Section 
5.0. The material of the pole piece was made of high-permeability iron and 
in fact, the pole pieces are the modified versions of the core pieces used 
for the study described in Section 3.0 and also for earlier work reported in 
Reference 4. As described in Reference 4, the teeth were .020 inch wide and 
separated by a gap of .040 inch and were cut to precision by electrical dis¬ 
charge machining. 

6.2.2 Test Rig 

The same test rig as described in Sections 2.0 and 4.0 was used for the 
present work. The modified core piece with the permanent magnet, shown in 
Figure 6-2, was fixed to the top plate, and the lower core piece (that has 
been used previously) with matching teeth was mounted on the bottom plate as 
shown in Figure 6-3. The energizing coils, shown mounted on the lower core 
piece, were not used in the permanent-magnet experiments described in this 
work. The gap and the transverse displacement were varied by moving the top 
and bottom plates respectively and were monitored by dial indicators and also 
by DC-LVDT'S. 
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The alignment of the core pieces was very critical in these tests because of 
the very small width of the teeth. Also, extreme care was taken to obtain 
the same air gap at both the legs of the core pieces. Since the transverse 
displacement was limited to 0.030 inch (one-half the pitch of the teeth), the 
motion was obtained by directly manually moving the bottom plate mounted on 
ball bushings with the help of a screw rod. A view of the test rig with the 
core pieces mounted is shown in Figure 6-4. 

6.2.3 Force Measurement 

The previously used, multi-directional quartz piezoelectric force transducers 
were used in the present study. The outputs of the force transducers were 
connected to two charge amplifiers through summing junctions. The outputs of 
the charge amplifiers yielded directly the normal and the transverse forces. 

6.2.4 Flux Measurement 

Flux measurement was performed by using a 10-turn search coil, mounted around 
the teeth of the lower core piece as shown in Figure 6-1. The search coil was 
designed to be removable from the magnetic circuit for the measurement of the 
absolute value of the flux. Care was taken to wind the search coil closely 
around the teeth to avoid any possible linking of the leakage fluxes with the 
search coil. 

The output of the search coil was connected to an integrating-type fluxmeter 
(LDJ Electronics Model 101). The fluxmeter output directly gave the number 
of lines of flux linking the 10 turns of the search coil. For purposes of 
obtaining plots, the fluxmeter output was directly connected to an X-Y plotter. 

In order to measure the flux at any required gap, the search coil is mounted 
around the magnet, and the two magnets are brought together until separated 
by the required air gap. The fluxmeter is reset to zero, and the flux-versus- 
gap curve is plotted as the two magnets are withdrawn from each other. To 


measure the absolute value of the flux, it is necessary to take the search 
coil away from the magnet to a distant point. 










With the permanent magnets in the magnetic circuit, the degaussing of the 
high-permeability material present in the circuit is not possible. In order 
to obtain repeatable values of fluxes and forces, the material of the core 
pieces was cycled through the flux changes before each experiment by bringing 
the upper core piece in close proximity of the lower core piece (gap of .005 
inch) and retracting it and then repeating this procedure about ten times. 

The increasing and decreasing flux curves did not deviate by any substantial 
amount after the above procedure was followed. 

6.3 Results and Discussion 

The experimental results will be presented in two sections, each section for 
a different magnetization level. A comparision of the experimental results 
with theoretical predictions is made in a separate section. 

6.3.1 One Magnet in the Circuit 

This section presents measurements of flux and forces as functions of gap or 
displacement when only one magnet (2 in. * 2-1/8 in. * .4937 in.) is used to 
energize the circuit. 

6.3.1.1 Flux Versus Gap . Figure 6-5 shows a plot of flux versus gap. The 
flux was measured with search coil in the position shown in Figure 6-1. Both 
increasing and decreasing flux curves were recorded. The deviation obtained 
between the two curves was considered negligible for the present work. The 
absolute value of the flux was obtained by moving the search coil from the 
magnetic circuit to a distant point. The flux reduces with increasing gap, 
but it appears that the decrease in flux is more gradual than found in earlier 
measurements when the circuit was energized electrically. 

6.3.1.2 Z-Force Versus Gap . The variation of Z-Force (attractive) versus 
gap is shown in Figure 6-6. The minimum gap reached was .005 inch. Fairly 
high normal forces are experienced at small values of gap. Here again, 
negligible deviation was noticed between force-gap curves for increasing and 
decreasing values of flux. 



































6.3.1.3 Flux Versus Displacement . The flux linkages through the magnetic 
circuit vary with transverse displacement because of the change in magnetic 
circuit reluctance with displacement. 

Figure 6-7 shows the changes in flux linkages at gaps of .005 inch and .010 
inch as the core pieces are displaced from the aligned position. The drop 
in flux linkages at large displacements appears to reduce significantly with 
increase in gap. Although measurement of normal forces was not made in the 
misaligned position, one can expect a similar drop in normal forces at large 


6.3.1.4 X-Force Versus Displacement . The transverse forces, which are 
stabilizing in nature, are plotted versus transverse displacement for variou 
values of gap ranging from .005 inch to .060 inch in Figure 6-8. The total 
transverse displacement traversed was .030 inch, corresponding to one-half 
the pitch of the teeth. The zero displacement in Figure 6-8 corresponds to 
a perfectly aligned position of the teeth. 


As the core pieces are transversely displaced, the X-Force increases indicat 
ing stability in the transverse direction, the X-Force reaches a maximum and 
then decreases with further increase in displacement. The maximum positive 
stiffness occurs at zero displacement. The transverse forces reduce in mag¬ 
nitude at increased gaps. It can, however, be noted that the flux linkages 
also decrease at increasing values of gap. A comparison between these exper 
mental results and theoretical predictions will be made in a later section. 


6.3.2 Two Magnets in the Circuit 

This section contains results for the case when two magnets of the same 
size as used before (and magnetized in the same direction) are used to ener¬ 
gize the circuit. Figures 6-9, 6-10, and 6-11 show curves of flux versus 
gap, Z-Force versus gap, and X-Forces versus displacement for the higher 
level of magnetization realized through two magnets energizing the circuit; 
these curves are very similar to those shown in Figures 6-5, 6-6, and 6-8. 
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The values of fluxes and forces for the two-magnets case do not appear to 
differ much from the corresponding values for the one-magnet case. At a 
measured flux level of 168,000 lines for the one-magnet case, the actual flux 
density in the teeth (considering the cross-sectional area of the teeth only) 

2 

is 1.84 Wb/m . The material of the pole pieces in the neighborhood of the 
teeth is, therefore, very nearly saturated. Under these circumstances, an 
addition of one more magnet to double the ampere-turns in the circuit does 
not bring about any noticeable change in the flux level. Since the flux 
level does not change for the two-magnets case, both the normal and the trans¬ 
verse forces remain approximately the same. 

6.3.3 Comparison with Theoretical Predictions 

The results reported above will be compared with the analytical prediction of 
forces described in Section 2.0. 

6.3.3.1 Force Calculation and Comparison . Table 6-1 shows the values of 
calculated normal and transverse forces corresponding to the experimental 
results described in Section 3.1. The flux variation with respect to gap was 
known from experimental results described in Section 3.1.1. Knowing the flux, 
the normal and transverse forces were calculated using the theoretical results 
described in Section 2.0. The numerical values of various quantities calcu¬ 
lated are shown in Table 6-1. The change in air-gap permeance and the corres¬ 
ponding change in mmf with lateral displacement of the teeth have been 
considered. Also, for cases when the flux changes were measured with lateral 
displacement, the corresponding changes in mmf have been taken into account. 
Gaps larger than .020 inch = 1.0) have not been reported in Table 6-1 be¬ 
cause of the small values of transverse forces and the limitations placed by 
the theory (-^ > .5). 

The above theoretical calculations are compared with experimental results in 
Figures 6-6 and 6-12 through 6-15. An excellent agreement between theoretical 
predictions and experimental results can be seen in Figure 6-6, which shows 
the variation of normal force with gap. 
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Figures 6-12 through 6-15 show measured values of transverse force at various 
values of gap versus transverse displacement compared with theoretical pre¬ 
dictions (marked X). The range of expected validity of the theory > .5) 

has been clearly marked in these figures. Excellent agreement between theory 
and experiments can be seen in the above figure within the range of validity 
shown. At large values of transverse displacement, the theory deviates from 
experiments significantly. Keeping in mind the assumption of uniform flux 
under the teeth made in the theory, this deviation is not surprising. Also, 
since in practical bearings the lateral displacements are small, the above 
fact does not severely limit the applicability of the theory. 

6.3.3.2 Flux Calculation . In the previous section, forces were calculated 
based on measured values of flux, and good agreement between experiments and 
theory for predicting forces was witnessed. 

Attempts were made using elementary analogies for the magnetic circuit to 
predict flux as a function of gap. The simple computations of the reluctances 
of the material of the pole pieces, the air gap, and the magnet (as a function 
of the gap and the flux level) were inadequate in predicting the flux-gap 
characteristics. It was concluded that, without accounting for the leakage 
fluxes, accurate prediction of flux-gap characteristics is not possible. 
However, leakage fluxes can be computed using existing, sophisticated finite- 
difference schemes. It was not the intent of the present research to include 
such a computation. 






7.0 CONCLUSIONS 

Subsequent to the magnetic bearing research described in an earlier report [4], 
several aspects of both closed and open magnetic circuits were studied, and 
these have been reported here. 

The theory that was developed for predicting forces in closed-circuit-type 
magnetic suspension and reported earlier [4] had been validated under static 
conditions only. Dynamic measurements of forces and fluxes were made for the 
same magnetic circuit to determine any effects of transverse vibration. 

Although the measurements discovered some degradation of transverse forces 
with frequency, the phase angle measured between the force and the displace¬ 
ment indicated very little damping in the system. Based on these findings, 
it appears that any damping required in such magnetic suspensions may have to 
be provided by external means only. 

Towards the study of open-circuit-type magnetic suspensions, a theory was de¬ 
veloped for predicting repulsive forces (normal and transverse) between per¬ 
manent magnets. This theory was based on the mathematical modeling of permanent 
magnets by their equivalent current densities. The repulsive force between 
two rectangular parallelopiped-shaped magnets were also measured. An excellent 
agreement was found between measured forces and theoretical predictions. 

Furthermore, the theory that was developed for the analysis of closed-circuit 
suspensions referred to earlier was verified experimentally for the case when 
permanent magnets were used as a battery to energize the magnetic circuit. 

Based on the research performed thus far, validated theories for predicting 
forces at the magnetic interface for both open- and closed-circuit systems 
appear to be at hand. However, the validation has been performed only for 
plane geometry. It is therefore necessary to validate the above theories for 
a rotating system with circular teeth. 
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